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We present the optimal local protocol to distill a 
Greenberger-Horne-Zeilinger (GHZ) state from a single copy 
of any pure state of three qubits. 
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Because of its relevance in quantum information the- 
ory, entanglement has been attracting considerable atten- 
tion in recent years. The initial efforts, mainly devoted 
to acquire a quantitative description of bipartite quan- 
tum correlations, led to the identification of the measures 
governing pure-state entanglement both in the asymp- 
totic and non-asymptotic (2) regimes. As a result, 
it is nowadays known how a bipartite system prepared 
in a pure state can be optimally manipulated under the 
restriction that only local operations on the subsystems 
aided with classical communication (LOCC) are allowed. 

Once a relatively complete command on pure-state en- 
tanglement had been achieved for bipartite systems, ef- 
forts have moved quite recently to address the somehow 
more intrincate case of tripartite entanglement, for which 
a system consisting of three two-level subsystems — that 
is, a three-qubit system — is the simplest scenario. 

There are two distinct ways in which three qubits can 
be entangled ||, in the following sense. Let us iden- 
tify as equivalent all pure states that can be reversibly 
interconverted, with some finite probability of success, 
when the parties are allowed to perform LOCC. Then, 
whereas all entangled pure states of two qubits are equiv- 
alent to the Einstein-Podolsky-Rosen (EPR) state Q, 
l/\/2(|00) + 1 1 1) ) , it turns out that truly tripartite pure- 
state entanglement of three qubits is either equivalent to 
the GHZ state B, 



|Gffz> = -^(|ooo> + |iii», 



(i) 



or else to the W state |), I/V3(|00I) + |010) + |100)), 
with these two states being completely inaccessible from 
each other by means of LOCC. 

Although it is genuinely tripartite, the entanglement of 
the W state only maximizes two-qubit quantum correla- 
tions fa] . As far as three-qubit correlations are concerned, 
the GHZ state appears as the maximally entangled state 
B: it violates Bell inequalities maximally and it maxi- 
mizes the mutual information of local measurements, ft 
is also the only state from which an EPR state between 
any two chosen qubits can be obtained with certainty. 



In this Letter we are concerned with the distillation of 
a GHZ state from an arbitrary pure state of three 
qubits (see also 0). We present the optimal strategy to 
distill a GHZ state starting from a single copy of \ip). 
That is, we provide the local protocol that allows the 
three parties to transform the state 1^) into a GHZ state 
with the greatest a priori probability of success compati- 
ble with LOCC. Marginally, we also discuss the best ap- 
proximate transformation g of |^} into \GHZ) by means 
of LOCC, which turns out to simply consist in a concate- 
nation of local unitary transformations. Obviously, these 
results may be of practical interest in any situation where 
three parties wishing to have a GHZ state only share the 
state \ip), and where for some reason — for instance be- 
cause they are in sufficiently distant locations — the re- 
striction to perform only LOCC is not merely academical. 
From a more general perspective, several difficulties have 
been faced which stem from characteristic features of en- 
tanglement involving more than two subsystems, such 
as the non-existence of the Schmidt decomposition (that 
is, from the fact that not all three-qubit states can be 
brought, by means of local unitary transformations, into 
the form a|000) + /3|If I), see ||). In this sense, our work 
explores the nature of entanglement beyond the bipartite 
case, while providing useful tools for its study. 

Let us first briefly review how, in the two-party setting, 
Alice and Bob can optimally distill a generic pure state 
of two qubits into an EPR state. The initial state can be 
transformed, by means of local unitary operations on Al- 
ice and Bob's sides, into a reference state whose Schmidt 
decomposition is 

a|00)+j8|ll) > (2) 

where a > (3 > and a 2 +[3 2 = 1. Then, following M, the 
diagonal operator M = — |0X0| -I- |1}(1| — corresponding 
to a POVM on, say, Alice's part — transforms state 
(||) into an EPR state with probability 2f3 2 . Finally, that 
such distillation protocol is optimal [flT| can be inter- 
preted in terms of the entanglement monotone E% = 1 
f§. Indeed, if a protocol would distill the EPR state 
with greater probability, this would contradict the non- 
increasing character of E2 under LOCC. 

Similar steps will be followed in order to find which 
is the optimal distillation protocol of a GHZ state in a 
three-qubit system: (i) the Schmidt decomposition will 
be replaced with a convenient two-term product decom- 
position; (ii) a subclass of distillation protocols, namely 
those in which each party is only allowed to perform one 
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positive operator-valued measurement (POVM) and only 
one of the overall outcomes is a GHZ state, will be dis- 
cussed and optimized; (Hi) the optimal of such one suc- 
cesful branch protocols (OSBP) will be proved to be the 
best distillation protocol by showing that the associated 
probability of success is an entanglement monotone p2| . 

Recall that a necessary and sufficient condition for 
the distillability of a GHZ from a three-qubit pure 
state \ip) is that the three reduced density matrices 
Pa =Trsc|' ! /')( ? /'|j Pb and pc have rank 2, and that the 
range of psc = r ^ r A\' l P){' l P\ contains two product vectors, 
say |&iCi) and I&2C2) |§- Only in this case \tp) admits a 
(unique) two-term product decomposition, 

\ip) = /ii|ai&ici) + p 2 e lip \a 2 b 2 c 2 ), fii > p 2 > 0, (3) 

where ip (z [0, 2ir) and the normalized local vectors 
l K 2) € C 2 satisfy 1 > (ki\k 2 ) > 0, k = a,b,c. In 
its turn, a distillation protocol consists of a scries of lo- 
cal POVMs, where the specific POVM to be performed 
at any time may be choosen depending on previous out- 
comes, the whole protocol having a tree structure. For 
a given branch of outcomes the initial state \tp) is trans- 
formed into a final state 



A®B®C\il>), 



(4) 



where the local operator A collects contributions from all 
POVMs performed by Alice in that branch, and similarly 
for operators B and C . Characterizing operator A by 
how it transforms \a\) and |a 2 ), 



A = a,\\a'i)(ai 



■ a 2 e ^ a \a 2 ){a 2 \ 



(5) 



where ai, a 2 , (a'^a^^O and (02)} is the biorthonor- 

mal basis to {|ai), |<Z2)} ; i-e. (o-i\aj) — <Jy, we ob- 
serve that the modifications introduced by A in \ip), 
A® 1 B $ lc|V) = ai/nK&iCi) + a 2 H2e i & + ^\a' 2 b 2 C2), 
concern the weights p\ and p 2 , the relative phase e %{p 
and the local scalar product (&i|a 2 ), whereas the other 
two local scalar products, (&i and (ci\c 2 ), remain un- 
changed. We can therefore readily conclude, from the 
uniqueness of the product decomposition (||) and the 
previous observation, that the only way for the result- 
ing state (0) to be proportional to \GHZ) is that (i) 
each local operator transforms the corresponding couple 
of non-orthogonal, local states into an orthonormal pair, 
|«2)} — > {|0), |1)}, k = a,b,c, that is, 



A = a 1 \0)(h 1 \+a 2 e^\l)(a 2 \., 
B=A|0X&i|+A 2 e i ^| 1 X^|, 

C = 7 1 |0)( Cl |+72e^|l)<c 2 |, 



(G) 



and that (ii) their combined effect modifies the weights 
Pi and p 2 to be equal, pi — > ouPc^iPi — \Jp/2, i = 
1,2, PI, where 



P ; 



(il>\A*A ® B*B <g> Ct(7|V) = 2(ai/3i7i^i) 



(7) 



is the probability that the distillation protocol follows 
that particular, successful branch. 

We have thus seen that the three parties must act on 
the system in order to distill a GHZ, because each of 
them must orthonormalize its local couple of states. The 
other relevant feature of the distillation process, namely 
making the relative weights equal, may be distributed in 
several ways among the parties. 

Our OSBP for distillation consists in each of the parties 
performing a unique, two-outcome POVM, and it is built 
in such a way that after each POVM one of the two pos- 
sible resulting states contains no three-partite entangle- 
ment at all, so that only one branch of the whole protocol 
succeeds in distilling a GHZ (see figure 1). In mathemat- 
ical terms this implies that if e.g. A in (g) is the success- 
ful operator in Alice's POVM, then the other operator of 
the local POVM, A, which satisfies A^A + A^A = 1 A , has 
rank equal to 1, so that the resulting state A^>1b^>^-c\' i P) 
is product in Alice's subsystem. Expressing the identity 
operator as 



1.4 



i,0 



(ai\aj)\ai)(aj\ 



and requiring that detfl^ — A' A] = 0, we find 
(l-a\){l-a 2 ) = {a 1 \a 2 ) 2 . 



(8) 



(9) 



Operators B and C are also accordingly constrained by 



(1-/3?)(1-/3|) = <6i|6 2 > 2 , 
(l-7i 2 )(l- 72 2 )Hci| C2 ) 2 . 



(10) 



To completely characterize our protocol, we will further 
impose that it succeeds with the greatest possible prob- 
ability over all OSBPs. The optimal OSBP probability, 



P(tp) = max 2(«i/3i7i^i) 2 , 



(11) 



— where the constraints aiftijipi = a 2 f3 2 r y 2 p 2 and (^|- 
|l(]) hold — reads 



p max 

r ~x>0 



1-1 1- 



4(1- ( ai |a 2 ) 2 ) 



/ 2 (*) 



1-1/1- 



4^|(1 -(6 1 |6 2 ) 2 )(l-<c 1 | C2 ) 2 ) 



/ 2 2 (*) 



(12) 



where 



h{x)= (x 2 + l)/x, 

f 2 {x) = (p\x 2 + 2pip 2 (bi\b 2 )(a\c 2 )x + p\)/x. (13) 

The maximization in ( ]l2] ) involves a polynomial equation 
of degree 6, which in general requires numerical calcula- 
tions. For illustrative purposes and later reference, we 
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will next consider two particular situations which can be 
solved analytically. 

Suppose, first, that the parties share the state 

=/ii|Q0ci) +/i 2 |llc 2 ), /ii>M2>0 (14) 

that is, with the local states already orthogonal in Alice's 
and Bob's subsystems. Then the previous optimization 
over all OSBP leads to a maximal probability 



p(0 1 ) = i- V 1 - iSiCi -(ci|c 2 ) 2 ). (15) 

This probability is two times the smallest eigenvalue 
of pc, which has a decreasing behaviour under LOCC. 
Monotonicity of this eigenvalue implies therefore that 
P(^) 1 ) is also the maximal probability for any general 
distillation protocol. It turns out that Claire — i.e. the 
party that has to orthonormalize its local states — is the 
only one that needs to perform a local POVM. 

The second particular case concerns states of the form 

\4> 2 ) = //i|06ici> +Ai 2 |16 2 c 2 ), Mi>M2>0 (16) 

that is with orthogonal states in Alice's subsystem. The 
maximum probability for OSBP is 

P{4> 2 ) = (l + 2 Ml /X 2 <& lC i|b2C2» 

J, 



/ 4/i^(i _ (& 1 |b 2 ) 2 )(l - ( Cl \c 2 ) 2 ~) \ 
V (l + 2p im (b lCl \b 2 c 2 )) 2 r [ ' 



Again, only Bob and Claire need to act on the system, 
with the corresponding POVMs satisfying 



Pi _ Hi jUi(&i|& 2 ) + A*2<ci|c 2 ; 7i 

7? 



/3f Mi 



P 2 /^(Mfo) + Ml( c l| c 2) 



(18) 



So far we have analyzed the optimal distillation pro- 
tocols under the constraint that only one branch leads 
to a GHZ. In what follows we will show that no distilla- 
tion protocol can succeed with probability greater than 
that for OSBP, P(ip). In order to do so, we will study 
the behavior of P{4>) under LOCC, to conclude that it 
is a decreasing entanglement monotone. That is, given 
the state \ip) and a sequence of local quantum operations 
that transform it into \ipi) with probability pi, we will 
show that 



(19) 



which means that the average probability to obtain a 
GHZ state from using several branches is not greater 
than when using just one branch. 

Although the set of transformations LOCC is very 
large and one should in principle check ( [l9| ) for any local 
protocol, we can use the fact that any such protocol de- 
composes into individual POVMs. Indeed, we need only 



prove the monotonic character of P(ip) under the most 

general local POVM on each subsystem. But, as a matter 

of fact, due to the symmetry of the problem it suffices to 

consider a general local POVM performed by one of the 

parties, say Alice. Notice, furthermore, that any POVM 

can be decomposed into a sequence of two-outcome 

POVMs Q. Let us then consider a two-outcome POVM 

with operators {N 1: N 2 }, nINi+N^N 2 = 1 A , applied by 

Alice. With some probability pi the resulting state will 
1 1 2 

be \ipi) = p i 1 Ni®\s® \c\i>)i an d then the parties 
can apply the optimal OSBP, with Ai® E>i® Ci being 
the corresponding successful branch (see figure 2(i)). We 
want to show that 



PW>PiP(^i)+p 2 P(^ 2 ). 



(20) 



Notice that the global action on Alice's side in or- 
der to distill a GHZ state can be reproduced by 
means of a single four-outcome POVM, namely by 
{AiNi, AiNi, A 2 N 2l A 2 N 2 }, where Ai are the operators 
satisfying A\Ai + A\Ai — 1a- The second and fourth 
operators disentangle the state. They are irrelevant for 
distillation, and we can join them together into some 
other operator R (figure 2(ii)). Both operators A\N\ 
and A 2 N 2 , being the last transformation Alice applies to 
its subsystem, have to leave the local states orthonormal 
and therefore can be written as 



AiNi = a M |0)(<fi| + a 2 , 4 |l)(a 2 |. 



(21) 



Consequently, this three-outcome POVM can be replaced 
with a two-outcome POVM with operators R and 



Q = y/<4 tl + «? i2 |0><ai| + y/c4 tl + «| )2 |l)(a 2 |, (22) 

followed — if the outcome corresponds to operator Q — 
by a diagonal, two-outcome POVM (figure 2(iii)), 



A = 



= |0)(0| + 



a 2 ■ 



Q 



1.2 



= |1)(1|. (23) 



a 



2.1 



*2 : 2 



Let \(f> 2 ) in (|lrj ) be the normalized vector after applying 
Q to Then we can translate condition (|2tj| ) into 



P(^ 2 )>9iP(0i) + 9 2 P(^), 



(24) 



where qi is the probability of the outcome related to Di 
starting from \4> 2 ), and \<p 2 ) is the corresponding final 
state. 

In order to finally check that (|24|) holds, we further 
notice that the action of any two-outcome POVM on a 
given state \ip) can be reproduced by a conditional se- 
ries of balanced, two-outcome POVMs, namely POVMs 
such that, for each of them, the probability of the two 
outcomes is exactly 1/2 |Tq ]. That is, we only need to 
analyze a diagonal POVM where the square of the diag- 
onal elements are x/ (2p\) and (1 — x)/{2p\) for the first 
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operator D\, where x € [2[if — 1, 2uf|, and their com- 
pletion to 1 for the second. Then ( p4| ) is an inequality 
for x that exhaustive numerical calculations have shown 
to saturate only for x = /xf, that is, when Alice acts 
trivially on her subsystem with Di = 1/^/21 a- Thus we 
conclude that the optimal probability for OSBP cannot 
be improved by any distillation protocol. 

Finally, we have also considered the optimal local ap- 
proximate transformation Q of \tp) into a GHZ state. 
That is, given the state \ip) and an arbitrary local pro- 
tocol which transforms it into the (possibly mixed) state 
pi with probability pi, we have looked for the maximal 
averaged fidelity between the output states pi and a GHZ 
state, 



F opt = mzxJ2Pi{GHZ\ Pi \GHZ). 



(25) 



As we argued in the bipartite setting g], we only need 
to consider final pure states {4>i,Pi} due to the mani- 
fest linearity of ( |25| ) in pi and the possibility of purifying 
any protocol that produces mixed states. Exhaustive nu- 
merical search shows then that the optimal approximate 
transformation consists in a deterministic transformation 
of \ip) into \ip') by means of local unitary transformations, 
exactly as in the bipartite case. 

We have addressed the optimal distillation of a GHZ 
state starting from a pure state of a three-qubit system. 
We have shown that, contrary to what happens in the 
bipartite case, tripartite distillation requires that each 
party performs at least one local POVM. This feature 
can be associated to the lack of Schmidt decomposition 
for tripartite pure states. An alternative decomposition, 
based on non-orthogonal product states, has proven very 
useful to indicate which modifications must be intro- 
duced locally in a pure state in order to distill it into 
a GHZ state, and to optimize such distillation. 
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A third condition to be satisfied is that ifi-\-ifi a +(pb-\-<Pc = 

mod 2n, although this can be accomplished by a simple 

local unitary transformation once the local states 

have been orthonormalized, and it does not affect the 

probability of a successful distillation. 

For instance, a three outcome POVM with operators 

{i?i}f =1 , yV E]Ej = 1, can be decomposed into a first 

two-outcome POVM with operators 



{ Ex, R =\ E ] 2 E 2 + E ] 3 E Z } 



(26) 



and a second two-outcome POVM, to be applied only if 
the outcome of the first POVM corresponds to the oper- 
ator 7?, with operators {E 2 R~ 1 , -E3-R -1 }. 
Consider an unbalanced, two-outcome POVM {Ei,E 2 }, 
eIE! + eIe 2 = 1, applied to the state with q > 1/2 
[1 — q] being the probability of obtaining the outcome 
associated with E x [E 2 \. Then the POVM 



{ 



£1, R 



^-1eIe 1 + eIe 2 } (27) 



transforms the initial state according to Ei with prob- 
ability 1/2. If the outcome corresponds to operator 7?, 
then a new 2-outcome, balanced POVM can be applied, 
which with probability 1 /2 will give one of the originally 
wished outcomes and with probability 1/2 will require a 
next POVM. The whole series of balanced POVMs can 
be easily seen to converge towards the original unbalance 
POVM. 
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FIG. 1. Scheme of a one succesful branch protocol (OSBP) 
for distillation of a GHZ state. 
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FIG. 2. As far as distillation of a GHZ is concerned, these 
three series of conditional POVMs are equivalent (see equa- 
tions (f|]|24j)) 



5 



